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ABSTRACT 

We describe and execute a novel approach to observationally estimate the lifetimes of 
giant molecular clouds (GMCs). We focus on the cloud population in the zone between 
the two main spiral arms in M51, i.e. the inter-arm region, where cloud destruction 
via shear and star formation feedback dominates over cloud formation processes. By 
monitoring the change in CMC number densities and ensemble properties from one side 
of the inter-arm to the other, we estimate the cloud lifetime as a fraction of the inter-arm 
travel time. We find that cloud lifetimes in M5I’s inter-arm are finite and short, i.e. 20 to 
30 Myr. Such short lifetimes suggest that cloud evolution is influenced by environment, 
in which processes are sufficient to disrupt GMCs after a few free-fall times. Over most 
of the region under investigation shear appears to regulate cloud lifetimes. As the shear 
timescale increases with galactocentric radius, we expect cloud destruction to switch 
primarily to star formation feedback at larger radii. We identify a transition from shear- 
dominated to star formation feedback-dominated cloud disruption through a change in 
the behavior of the inter-arm GMC number density. The signature suggests that shear 
is more efficient at completely dispersing clouds, whereas star formation feedback tends 
to transform the cloud population, e.g. by fragmenting high mass clouds into lower mass 
pieces. Compared to the characteristic timescale for molecular hydrogen in MSI, our 
short cloud lifetime measurements suggest that gas can remain molecular while clouds 
disperse and reassemble. We propose that galaxy dynamics regulates the cycling of 
molecular material from diffuse to bound - and ultimately star-forming - objects, and 
hence contributes to long observed molecular gas depletion times in normal disk galaxies. 

We also speculate that, in more extreme environments such as elliptical galaxies and 
concentrated galaxy centers, star formation can be suppressed when the shear timescale 
becomes so short that some clouds never have the opportunity to collapse and form stars. 


1. INTRODUCTION 

The lifetimes of giant molecular clouds (GMCs) 
set a natural limit to the timescale over which gas 
is converted into stars and thus potentially impose 
powerful constraints on simulations of galaxy for¬ 
mation and evolution. Yet very little is directly 
observationally known about cloud longevity. Re¬ 
cently, a wealth of observations at high spatial 
resolution and sensitivity in nearby galaxies have 
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stimulated renewed debate on the issue of cloud 
lifetimes, their relation to the molecular content 
of galaxies and the process of star formation at all 
spatial scales. 

The very fact that we observe clouds in spiral 
galaxies populating the area between spiral arms 
(in the so-called inter-arm region), has been taken 
as evidence that they are long-lived, surviving at 
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M33. Such short lifetimes are consistent with re¬ 
cent models of GMC formation and evolution that 


predict life t imes nearer 10-20 Myr (e.g. Dobbs 
et aU[20TTj [Dobbs et al. |201^ [Dobbs fc Pringle 
2013p. In Emulations with realistic distributions 
ot gas and stars, ISM heating/cooling and star for¬ 
mation, the combination of shear and star forma- 
tion feedback lead t o fairly rapid cloud disp ersal 
(Dobbs et al. 2012 Dobbs & Pringle 2013). In 
support of this picture, two recent studies of the 
cloud population in M51 conclude that clouds are 
far from ‘standard’ and do not obey the scaling 
relations normally taken to imply that they are 
virialized ob jects (Hughes et al. 2013b Colombo 
et al. 2014a). This suggests that clouds are not 
necessarily stable, long-lived entities. 

Whether clouds are long- or short-lived has 
strong implications for the efficiency of star for- 
mation observed across a range of spatial scales 
( Krum holz & Tan 2007 Feldmann & Gnedin 
201 1|, impacting our view of the balance between 
cloud-scale physical processes thought to regulate 
star formation, e.g. turbulent driving, magnetic 
fields and feedback in the form of mechanical and 


thermal energy (jVazquez-Semadeni et al. 

2003 
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]). In order to distinguish between 


static views of GMC evolution, a greater num¬ 
ber of direct observational estimates for GMCs 
lifetimes in systems beyond the Local Group is 
urgently required. To date, the high resolution 
data required to establish an empirical picture of 
the GMC lifecycle has only been available for the 
Milky Way, the Magellanic Clouds and M33, and 
hence the diversity of galactic environments for 
which there are empirical estimates of GMC life¬ 
times is very limited. 

In this paper, we describe a technique for esti¬ 
mating GMC lifetimes in star-forming disk galax¬ 
ies that is independent of their association with 
young stellar phenomena. We exploit the large 
number statistics provided by the PdBI Arcsecond 
Whirlpool Surve y (PAWS; Schinnerer et al.[|2013 


Pety et ^[2013 ) to monitor changes in the CMC 


population across the inner disk of M51, which we 
then interpret using our detailed knowledge of gas 
dynamics and high mass star formation across the 
PAWS field. The premise is a simple one: a short 
(i.e. less than an orbital period) cloud lifetime 
should manifest itself as a decrease in the num¬ 
ber of clouds from one side of the inter-arm to the 
other. The ratio of the initial number of clouds 
to the number of clouds ‘lost’ during inter-arm 
passage constrains the ratio of the cloud lifetime 
to the time to traverse the inter-arm. Since the 
travel time in the inter-arm increases with galac- 
tocentric radius the fraction of ‘lost clouds’ may 
exhibit a radial dependence: fewer clouds are ex¬ 
pected to survive over the full passage through the 
inter-arm at larger galactocentric radius. The sig¬ 
nature may also depend on the dominant cloud 
disruption mechanism (shear, star formation feed¬ 


back) within a given zone of the galactic disk. 

This paper is organized as follows. In § [^ we 
present our method to estimate cloud lifetimes. 
As we are primarily interested in cloud destruc¬ 
tion processes and the limits they impose on cloud 
lifetimes, we focus on the cloud population in the 
inter-arm region of disk galaxies, i.e. the zone in 
between the main spiral arms. Gompared to the 
high density arm environment, fewer clouds are 
expected to form in the inter-arm and, moreover, 
the dynamics in this region of the disk are easier to 
characterize. In § [2.4[ we summarize the strengths 
and weaknesses of our method compared to other 
approaches for estimating cloud lifetimes and dis¬ 
cuss its general applicability. In § [^ we apply our 
framework to M51, where we use the pattern of 
star formation and shear in M51’s inter-arm zone 
to separately study the mechanisms that act to 
limit the lifetimes of clouds. We present our mea¬ 
surements of cloud lifetimes in MSl’s inner disk 
in § We discuss and interpret radial trends in 
cloud^lifetimes and cloud destruction mechanisms 
in §[^ Finally, in §|^we relate our findings to pre¬ 
vious cloud lifetime estimates, emphasizing their 
impact on our view of cloud-scale star formation. 

2. THE MEASUREMENT FRAMEWORK 

In this section we outline our method to es¬ 
timate cloud lifetimes, presenting first a heuris¬ 
tic and then a more quantitative description of 
our approach. Since we aim to provide an anal¬ 
ysis framework for application to real data, we 
briefly discuss the connection between the ideal¬ 
ized clouds in our model and GMCs as identified 
in CO observations. We refer to the PAWS CO(l- 
0) survey to illustrate these general considerations, 
which are common to all cloud-scale extragalactic 
CO surveys to which our method can be applied. 
Later in § [^ we apply our technique to the inter¬ 
arm cloud population in M5I to derive an estimate 
of the characteristic GMC lifetime. 


2.1. Cloud formation and destruction in disk 
galaxies 

The key element for measuring cloud lifetimes 
with our method is a survey of molecular clouds 
in the inter-arm zone between spiral arms. This al¬ 
lows us to track azimuthal variations in the inter¬ 
arm GMC number density, which we hypothesize 
reflect the balance between cloud destruction and 
formation processes within the inter-arm zone. 
This is motivated by our analysis of GMCs iden¬ 
tified in the PAWS survey of CO(1-0) emission in 
the central 9 kpc of M51 (Colombo et al. 2014a) 
which showed that the nuniber of inter-ami clouds 
in the inner disk of M51 decreases from the down¬ 
stream of one spiral arm (henceforth “zone I”, see 
Figure 1), to the upstream region of the other arm 
(“zone II”). The cloud mass spectrum also changes 
in shape across the inter-arm, consistent with the 
idea that the cloud population evolves across the 
inter-arm due to disruptive processes such as shear 
and feedback (Colombo et al. [[2014a). 
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Hereafter, we adopt the following simple picture 
for the passage of clouds through the inter-arm 
(shown in Figure[^: after leaving the downstream 
of one spiral arm, clouds enter zone I and proceed 
to zone II, which is located upstream of the next 
spiral arm. Under the assumption that clouds 
follow rou ghly circular paths through the inter¬ 
arm (see § |3.2[ ) the current azimuthal position of 
a cloud provides a measure of the time that has 
passed since it entered the inter-arm. Using the 
galaxy’s orbital period as a fiducial clock, we can 
thus connect eyolution in the number of clouds 
from zone I to zone II with an estimate of the 
characteristic GMC lifetime. 

As clouds pass from zone I to zone II in the inter¬ 
arm, we expect them to be susceptible to two pri¬ 
mary destructive mechanisms, shear and star for¬ 
mation feedback. Shear acts on all clouds in galac¬ 
tic disks undergoing differential rotation, and thus 
may be considered to set the ‘natural’ cloud life¬ 
time in the absence of other destructive processes. 
Star formation feedback should also contribute to 
cloud disruption, although the timescale associ¬ 
ated with destruction via feedback may depend on 
the specific energetics of the star formation event 
and a cloud’s proximity to the star-formation ac¬ 
tivity. 

These two processes may lead to different sig¬ 
natures in the variation of GMC number density 
with azimuth (although our estimate of the GMC 
lifetime does not rely on it). Shear, for exam¬ 
ple, may l ead to mass loss and/or c omplete cloud 
dispersal (|Dobbs & Pringle ||20I3|). In the pres¬ 
ence of shear, our basic expectation is that the 
number of clouds and their total combined mass 
smoothly decreases from zone I to zone II, as ma¬ 
terial in clouds is returned to the ISM. By con¬ 
trast, feedback, which can disrupt clouds as well 
as sweep up and co mpress interstellar gas (e.g. 


Dawson et al. 2013), is expected to be spatially 
localized and coincident with tracers of high mass 
star formation. In this case, mass loss may be 
less gentle and poten tially involve cloud splitting 
(Dobbs et al. 2011), thus forming new clouds. 
At the same time, feedback acting externally can 
produce converging flows , possibly merging pre¬ 
existing cloud fragments (Dobbs et al. 2012) and 
increasing cloud masses. As a result, we expect 
feedback-dominated zones to exhibit a transfor¬ 
mation of the cloud population (i.e. an exchange 
of mass between high- and low-mass clouds) near 
sites of star-formation, rather than a steady de¬ 
cline in both the number and combined mass of 
clouds. 

In either case, these considerations suggest 
that some modes of cloud destruction can be 
accompanied by an increase in clouds (usually 
an increase in the number of low-mass clouds), 
and not solely a decrease in the total number of 
clouds. In our model of GMG evolution in the 
inter-arm region, we therefore consider both sinks 
and sources, i.e. cloud destruction and cloud 
formation. 


arm zone I zone II arm 



time 


Fig. 1. — Polar coordinate representation of the inter¬ 
arm cloud population (blue ellipses) between two spiral 
arms (blue diagonal lines). The two solid black lines mark 
the boundaries of the inter-arm environment. The black 
dashed line shows the midpoint that divides the inter-arm 
into zone I and zone II. The distribution of star formation 
events (stars) qualitatively follows the observed pattern of 
star formation in M51 (see Figure®: orange stars repre¬ 
sent star formation located within me spiral arm environ¬ 
ment, while red stars indicate inter-arm star formation. In 
this representation, the trajectory of clouds through the 
inter-arm is in the horizontal direction (as indicated by the 
arrows). 


2.2. Measuring the cloud lifetime 

We now present a more formal description of 
our method for estimating the characteristic GMC 
lifetime. Although we cannot directly follow the 
evolution of clouds with only a single observational 
snapshot, we can adopt an Eulerian representa¬ 
tion of the gas flow given a sufficient number of 
clouds in our sample. This allows us to replace 
measurements of clouds throughout their lifetimes 
(Lagrangian representation) by statistical cloud 
measurements as a function of position in a spiral 
galaxy (in our case M51). In this way, for a cloud 
population of initial size Nq at time to that un¬ 
dergoes subsequent evolution, we can estimate the 
cloud lifetime t using a measurement of the cloud 
population at some later time U- Here we assume 
that a single lifetime r (the statistical average in 
the population) applies to all clouds, rather than 
a distribution of cloud lifetimes. In this case, if 
the population is reduced to A^i by time U, then 
the implied loss rate of — (A^o — A^i) / (^i — ) leaves 
no clouds left after time r, i.e. 


No- 


Nq-Ni 
ti — to 


T = 0. 


( 1 ) 


Although we assume that the population overall 
suffers losses, we let the rate —{No — A^i)/(ti — 
to) represent the effective rate of change in the 
population in the presence of sources and sinks, 
assuming that the loss and gain (or growth) rates 
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are independent of time, i.e. 


Nq-N, 
ti — to 


= {loss — gain) = Nq 


'^true 


• grow 


( 2 ) 

Here we specifically let Tgrow represent the time 
to increase the population by Nq at the current 
growth rate and take Ttme to be the time it would 
take to reduce the initial population at the given 
loss rate to zerOji.e. the true cloud lifetime. 

Eqs. Q and Q together imply 


1 


1 


'^true 


' grow 


( 3 ) 


and we see that when Tgrow»Ttrue-i the cloud life¬ 
time can be approximated by r. 

In principle, measurements Ni and Nq can be 
made in any two adjacent zones between which 
there is a well-dehned travel time ti — to. Since 
we use clouds sampled at different spatial loca¬ 
tions along a common trajectory as proxies for the 
evolution of a single cloud, the travel time between 
two zones separated by angular distance Acj) repre¬ 
sents the Lagrangian-equivalent time for a cloud 
to traverse that distance, from one zone to the 
next. 

To estimate the travel time requires knowledge 
of the angular velocity V^. In the inter-arm re¬ 
gion, where non-circular motions are negligible, 
clouds follow a circular path and can be ap¬ 
proximated^^ the circular velocity in the disk 
(see also § 3.3 for detailed discussion). Position¬ 
ing the measurement zones in this region, to avoid 
the spiral arms, leads to more reliable estimation 
of ttrav than when zones overlap with spiral arms, 
where non-circular motions are present and must 
be taken into account. Focussing on the inter-arm 
also offers a more direct measure of Urue, since 
clouds here are expected to undergo relatively less 
formation compared to, i.e. in the spiral arms. 

Given large enough numbers of clouds, in prin¬ 
ciple more than two zones in the inter-arm can 
be considered, providing more than one indepen¬ 
dent measure of the lifetime. But to maximize the 
number of clouds in each zone, we recommend a 
single lifetime measurement (per inter-arm) using 
two zones that together span the entire width of 
the inter-arm. In this case, eq. 0 simplihes to 


^travel 

2 Ni - Nil' 


( 4 ) 


Here, the cloud population spanning some (az¬ 
imuthal) area At/) at a snapshot in time is di¬ 
vided in half, into two populations, Nq=Nj and 
Ni=Nii. Measurements in the two halves (zones) 
are thus separated in time by the equivalent time 
in a Lagrangian representation for a test cloud to 
cross half the distance At/) (i.e. ttrai;e//2 where 
ttravei is the time to cross the full distance). 

In terms of the lost fraction Fiost = {Ni — 
Nii)lNj eq. 0 becomes 


ttravei 1 
2 Fiost ’ 


( 5 ) 


and we expect short cloud lifetimes to result in 
higher Fiost, or fewer clouds in zone H compared 
to zone I. 

In our recommended application of eq. 0 , mea¬ 
surements of the cloud number density in a spiral 
galaxy occur in the two halves of the inter-arm, 
zone I or zone H, so that ttravei is the length of 
time to cross the distance A^ between spiral arms 
at their present location[0 The only requirement 
is that the azimuthal distance At/) should be large 
enough that the corresponding travel time spans 
enough of a cloud’s evolution to be sensitive to 
factors that limit its lifetime. Note that the travel 
time defined in this way is shorter than the full 
length of time for a cloud to pass from one spiral 
arm to the next, tsp = 2-^/m{VL — Hp) where Hp is 
the so-called pattern speed of the spiral and m is 
the number of arms. 

In practice, relating evolution in the number of 
clouds from one side of the inter-arm (A/) to the 
other {Nil) with cloud lifetimes using eq. 0 
requires three key components: (1) a catalog of 
cloud positions and properties that is complete to 
a well-determined sensitivity limit; (2) an accu¬ 
rate measure of the rotation curve (to estimate 
the inter-arm travel time as a function of galac- 
tocentric radius); and (3) an estimate of radial 
streaming motions (to place a bound on the radial 
excursion of clouds in their orbits as they travel 
through the inter-arm). Each of these should be 
readily available for most (nearby) galaxies with 
existing state-of-the-art and future molecular gas 
surveys. In § [^ we provide an explicit example of 
how we estimate GMC lifetimes in M51 by apply¬ 
ing our method to the PAWS data. 


2.3. Clouds in observations and in theory 

The application of our method to actual obser¬ 
vational data involves an obvious but important 
assumption, i.e. that a galaxy’s GMC popula¬ 
tion is accurately described by the cloud catalog 
constructed from observations. In practice, this 
means that the cloud lifetime strictly applies only 
to objects that are well-represented by the ob¬ 
servational data, which depends on observational 
characteristics such as the resolution and sensitiv¬ 
ity of the survey. 

For example, the PAWS data that we analyze in 
§|^has a spatial resolution of ~ 40 pc and a spec- 
trm resolution of 5 km s“^, with an RMS noise 


of 0.4 K pe r channel (Pety et al. 2013 Schinnerer 
et al.||2013). Assuming a Galactic conversion fac¬ 
tor from Go luminosity to molecular gas mass, 
a 10^ Mq GMC in M51 should therefore be de¬ 
tected by the PAWS survey with 5a significance. 


Although eqs. 0 and 0 can be applied, in principle, 
to galaxies without ^rong spiral patterns (or even outside 
the inter-arm) provided an appropriate ttravei can be de¬ 
fined, in practice the estimated lifetime is highly uncertain: 
the cloud populations of fiocculent galaxies or those with 
weak, multi-armed spirals tend to be sparser than in galax¬ 
ies with strong, well-defined spiral arms, which provide a 
source of new clouds. Without large numbers of clouds, 
the method cannot be reliably applied. 
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This then translates into a completeness limit for 
the cloud catalog constructed from these (or sim¬ 
ilar) observations through decomposition of the 
CO emission, i.e. with CPROPS (see Colombo 
et al. 2014). Our definition of ‘cloud’ thus applies 
only to emission above the sensitivity limit when 
this emission is moreover connected in position- 
position-velocity (pp v) s pace (as assumed in the 
decomposition; see § 3.1). 


Although the lifetime we measure is arguably 
sensitive to the identihcation and de comp osition 
that dehnes clouds (as discussed in § 2.3, requir- 


2 . 3 . 1 . The cycling of molecular material 


ing, in particular, coherent structures in ppv space, 
it is independent of the actual cloud mass (and 
how it varies with time). We therefore avoid sev¬ 
eral of the uncertainties associated with mass esti¬ 
mation, from the definition of the cloud boundary 
itself to the conversion of the surface brightness 
of a given tracer to gas column density (e.g. the 
CO-to-H 2 conversion factor). 


Some part of the CO emission detected in galax- 


2 . 4 . 1 . Relation to star formation 


by PAWS ( 

Pety et al. 

2013 

). Roughly half of all 

The method is also independent of the presence 

CJU emission surveyed 

throughout the TAVVy held 

of star formation, which can often be difficult to 


of view is compo sed of clouds (dropping to ^40% 


in the inter-arm; Colombo et al. 2014a). The re¬ 
maining CO emission, which we henceforth refer 
to as part of the ‘non-cloud medium’, is partly 
resolved, but not organized into coherent ‘cloud¬ 
like’ en tities, and partly structured on much large r 
scales dColombo et ai||2014a| |Pety et al.||2013|). 
Regardless ot the nature ot this emission, when 


determine (due to strong dust extinction or sen¬ 
sitivity limitations). Our cloud lifetime measure¬ 
ment thus avoids such ambiguities and is neither 
limited to the onset of star formation nor sensitive 
to the duration of star formation. In principle, the 
measured lifetime can extend beyond the star for¬ 
mation event, as long as the cloud remains more 
massive than lO^M©. 


clouds evolve outside the definition of a cloud cat¬ 
alog we assume that they ‘die’: the measured life¬ 
time marks the time either when clouds are dis¬ 
persed back in to the non-cloud medium or when 
their masses and sizes fall below our detection 
threshold]^ Similarly, we assume that a cloud 
is not formed until it grows beyond the detection 
threshold of the dataset. 

2 . 3 . 2 . Phase changes and continuity 

In general, clouds could also ‘disappear’ as a 
result of phase changes. Whether clouds cata¬ 
loged with a particular tracer represent a complete 
sample (and thus supply a useful tracer of cloud 
lifetimes in an Eulerian representation) depends 
on whether the gas phase being traced undergoes 
transitions to other ph ases (as considered in the 
case of M5I in section § |3.5| ). 

Where the ISM is molecule-dominated, as in the 
central area of M51 covered by the PAWS field 
of view, we can use CO as our primary tracer of 
the gas including its cloud entities. In other in¬ 
stances (in other galaxies with different balances 
in their atomic and molecular phases), variation 
in the mass in a given phase from one side of the 
inter-arm to the other may imply that additional 
gas phases must be accounted for, or that the CO- 
to-H 2 conversion factor varies azimuthally. Before 
applying our method, we therefore recommend 
first assessing mass continuity across the inter-arm 
as considered later in § |3.5[ Note, though, that 
continuity is not required by our metho d, ev en if 
it may provide a useful diagnostic; see § |2.4| 

2 . 4 . (Relative) Strengths and weaknesses of the 
method 


2 . 4 . 2 . Population Growth 

Our method is less powerful for strictly esti¬ 
mating cloud lifetimes when the cloud popula¬ 
tion is also characterized by growth due to, e.g. 
cloud formation or transformation, which may 
often accompany star formation (and associated 
feedback). When clouds are introduced into the 
population between zone I and zone II, the mea¬ 
sured Float will underestimate the true fraction of 
‘lost’ clouds. Eq. ([^ thus in general provides 
only an upper limit on the cloud lifetime, unless 
the number of clouds added to the population is 
far exceeded by the number of clouds lost (i.e. the 
timescale Tgrow to gain clouds is much longer than 
the true cloud lifetime Ttme)- 

Under most circumstances in which the (inter¬ 
arm) population contains growth, we expect the 
timescales for population growth and loss to be 
similar, i.e. clouds emerge from the destruction or 
transformation of existing clouds. When compara¬ 
ble numbers of clouds are gained and lost between 
zone I and zone II, there will be no evolution in 
cloud number. As Float approaches zero, r will 
deviate strongly from the true cloud lifetimej^ In 
this scenario, eqs. and (§ can provide an esti¬ 
mate for the population growth timescale, assum¬ 
ing a specific cloud lifetime. 

2 . 4 . 3 . Other assumptions of the method 

Our assumption that the rates of gain and loss 
in the population are independent of time is one 
of necessity, but it should not be entirely unre¬ 
alistic, e.g. when describing losses due to shear. 
(Changes in the galaxy’s mass distribution and 
hence orbital velocities will typically vary over 


We note that the total emission in the cube does not 
support the presence of significant numbers of clouds below 
10®A4q, extrapolated from the low-mass end of th e mass 
spectra of cataloged clouds ([Colombo et al.|2014a^. 


For populations predominantly undergoing growth 
(with comparatively few losses) Fiost<0, and eq. nat¬ 
urally yields a lower bound on the timescale with which 
clouds are formed. 
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much longer timescales than an inter-arm travel 
time.) For other mechanisms, both the rate itself 
and how quickly it changes could depend on cloud 
properties, such as mass or surface density; if the 
process of star formation is limited to clouds of 
a certain type, then destruction via feedback will 
only act on this subset of the population. With a 
large enough sample of clouds, it should be possi¬ 
ble to separate clouds into subpopulations in order 
to identify whether such a dependence might exist. 

While these factors can be easily integrated into 
alternative approaches, e.g. in which cloud life¬ 
times are estimated by fitting tailored models mix¬ 
ing cloud formation and destruction to cloud mass 
spectra, our simple (reductive) approach affords 
relative model-independence, requiring only mea¬ 
surements based on observables..pdf 


3. APPLICATION TO M51 

3.1. Cloud-decomposed inter-arm emission 

Our estimate for the cloud lifetime in MSl’s 
inter-arm relie s on the catalog of G MCs identi¬ 
fied in PAWS (Colombo et al.|j2014a|. This cata¬ 
log contains ^ 15UU clouds, and includes measure¬ 
ments of the cloud position, size, linewidth, peak 
brightness, integrated CO luminosity and dynam¬ 
ical mass. Along with th e spiral arms, molecula r 
ring and nuclear bar (see Colombo et al. 2014a), 
the inter-arm is one of four main dynamical envi- 
ronments in M51’s inner disk. Just over 500 clouds 
populate this region in the disk, defined as the area 
between the environment of the two main spiral 
armand extending from R=1.3 kpc, from the 
‘center environment’, to the edge of the field of 
view (see Figure]^. This large number of clouds 
{Nciouds > 500) is required to accurately monitor 
changes in the cloud number density as a function 
of azimuth and galactocentric radius within the 
inter-arm environment. 

Inter-arm clouds are split into two populations, 
one on either half of the inter-arm as defined by 
Colombo et al. (2014a). The inter-arm sits be¬ 
tween either of the two main spiral arms, whose 
widths are determined via observed gas kinemat¬ 
ics as is, therefore, the location of the inter-arm 
mid-point itself. The distinction between zone I 
and zone II is illustrated in Figure showing the 
polar-coordinate projection of clouds in the inter¬ 
arm to the south relative to the two spiral arms. 
Although the particular sorting of clouds into zone 
I and zone II can be sensitive to the exact location 
of the inter-arm midpoint, quantifying our uncer¬ 
tainty in this location serves as a way of evalu¬ 
ating the dominant uncertainty on cloud lifetimes 
measured with our technique (as defined in section 
4.1). As there are two main spiral arms in M51, 
there are are two inter-arm regions. Following the 
sorting of clouds into zone I and zone II within 
each inter-arm (into four populations), their num¬ 
ber densities are determined and then combined 
together (into two main populations). 


The spiral arms have a finite, kinem aticallv- 
determined angular width ([Colombo et al.|2014a|. 


According to Colombo et al. (2014aI, the bo- 
sensitivity of the PAWS dataset is 1.2 x 10^ Mq 
throughout most of the field of view, but an in¬ 
crease in the noise toward the edge of the field in¬ 
creases the completeness limit to 3.6x 10® Mq. We 
therefore choose to include all cataloged clouds, 
which extend down to 1.2 x 10® M©, in our anal¬ 
ysis at all radii except at i?>70”, where we keep 
only clouds more massive than 3.6 x 10® Mq given 
the change in noise pattern. We find very little 
change in our results with small changes in this 
threshold, as most clouds are well above this mass. 


3.2. Cloud trajectories 

According to our definition of the inter-arm en¬ 
vironments, clouds leave the spiral arm down¬ 
stream and enter zone I and then proceed to zone 
II located upstream of the next spiral arm. The 
trajectory of clouds during their inter-arm passage 
is expected to be roi^hly circular (e.g. straight 
left to right in Figure^, but it is important to ver¬ 
ify that clouds do not drift significantly in galac¬ 
tocentric radius. In the event of large radial ex¬ 
cursions, changes in cloud numbers from zone I to 
zone II could be explained entirely as clouds pass¬ 
ing from zone I in one radial bin to zone II in the 
neighboring bin, or even passing directly between 
radial bins within zone I itself. 

In general, clouds in spiral galaxies are subject 
to large non-circular motions as they orbit in the 
non-axisymmetric potential. But once a cloud ex¬ 
its the spiral arm, where these motions are largest, 
the cloud unde rgoes nearly cir cular motion in the 
inter-arm (i.e. Roberts 1969). We can estimate 
a cloud’s radial excursion about its circular orbit 
as roughly the size of the epicyclic radius, Vr/n, 
where k is the epicyclic frequency and Vr is the 
size of radial streaming motions. These motions 
are modest in the inter-arm of M5I compared to 
the spiral arms. But we can place a conserva¬ 
tive estimate on the radial distance a cloud will 
traverse during its inter-arm passage by adopting 
the maximum val ue of radial inotion s in the spiral 


et al. 

2014b 

together with the radially varying n 

obtained by 

Meidt et al. 

(2013). We find that the 


radial excursion is everywhere only ~ 4% of the 
orbit circumference (even at a radius i?=4 kpc, 
this is at most 300 pc). Clouds can thus be safely 
assumed to follow circular paths as they cross the 
inter-arm environment. 


3.3. Travel time in the inter-arm 

Given that clouds in the inter-arm follow cir¬ 
cular paths, the time it takes for a test cloud 
to travel the current distance spanned between 
any two of m spiral arms can be estimated as 

ttravel = {‘2TT - m9arTn)Rm~^V^^^, where Oarm IS 

the angular width of the spiral arm and Viarm 
is the azimuthal inter-arm velocity. In M5I 
to= 2, since there are two spiral arms, and we 
approximate Viarm^Vrot (he. little inter-arm 
streaming) so that ttravei=r'R/Vrot- This travel 
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Fig. 2.— Map of in ter-arm cloud positions (blue) extracted from the PAWS survey of CO(l-O) emission in M51 by 
[Colombo et al.| | |2014a| |. Clouds in zone I of the inter-arm are marked with filled circles, while zone II clouds are shown 
with open circles. Cr ay contours highlight t he position of the spiral arms traced in CO. Contours of Ho; emission across 
the PAWS FOV (fromjSchinnerer et al.|2013| assuming thejGutierrez et al.|||2011|l stellar continuum correction of the HST 
ACS image) are shown color-coded by environment: the zones downstream and upstream in the inter-arm (red) and the 
arm and center environments (yellow). Concentric black rings mark radii i?=40”, 60”, 80” and 100”. 


time is good to within ^ 4 Myr, accounting 
for the assumed finite, uniform spiral arm width 
'u;=i?^arm"^300 pc. Note that although we write 
torb=‘^'^Ryrot=‘^ttravei throughout, this Underesti¬ 
mates the true time to complete one (non-circular) 
orbit, which is 2{tarm+ttravel) with tarm=wV~^l^, 
where Varm represents transverse motions through 
the spiral arm pot ential that reach ^ 30 km s“^ 
(Meidt et al. 2013). As we are interested only in 
the time to travel the distance spanned by the two 
spiral arms (and its uncertainty) at each radius we 
simply adopt the circular vel ocity from th e rota- 
tion curve model derived by |Meidt et al. (2013) 
from the baryonic mass distribution in Mol. 


3.4. Radial variation in cloud 
formation/destruction 

M51 presents a unique opportunity to distin¬ 
guish between two main cloud dispersal processes, 
i.e. shear and feedback from massive star forma¬ 
tion. Over the radial range spanned by the PAWS 
field-of-view, our knowledge of the galaxy’s mass 
distribution, gas kinematics and global pattern of 
star formation indicates that there are two dis¬ 
tinct radial zones where we can expect the influ¬ 
ence of each of these processes to dominate. As 
shown in Figure most of the massive star forma¬ 
tion in M51 occurs along the spiral arms, limited 



Fig. 3. — Polar coordinate representation of Fig. high¬ 
lighting clouds in the inter-arm to the south. Clouds in 
zone I are marked with blue filled circles, while zone II 
clouds are shown with blue open circles. As in Fig. 0 gray 
contours highlight the position of the spiral arms traced in 
CO, while contours of Ho emission are shown color-coded 
by environment: the zones downstream and upstream in 
the inter-arm (red) and the arm and center environments 
(yellow). 


to radii i?>60”[3 and appears offset just down¬ 
stream of the spiral arms traced at high resolu¬ 
tion in the PAWS CO(l-O) map. The majority of 

The evident suppression of star formation along the 
inner spiral segmen t, between 40”<i?<60”, is discussed by 
[Meidt et al.j ( [2013| l. 
































Fig. 4.— (Top) Radial profile of inter-arm shear mea- 
sured by the bac kground Oort A estimated in M51 by 
|Meidt et al.|(|2013j. fBottoml Radial trends in 
(inverse gas depletion time Trfep)i measured by Meidt et| 
|al.| ||2013|l. Measurements are extracted in radia,! bins in 
wtiicn azimuthal averaging of Sgas (estimated from PAWS 
CO and THINGS HI) and Ssfa (measured from Ha and 
24 fim) runs from 0 to 27r (throughout the PAWS FOV; 
gray) or only across the inter-arm (black). Uncertainties 
are on the order of 30% l|Meidt et al.|2013f. We interpret 
regions at ij>60” where the star formation rate per unit 
gas mass is high as locations where feedback from massive 
star formation dominates cloud destruction. The zone in 
which the star formation is reduced but shear is relatively 
high marks the region where we expect cloud destruction 
primarly through shear. 


this star formation falls within our definition of 
the arm environment, but it persists into zone I 
of the inter-arm. Clouds at these galactocentric 
radii thus appear susceptible to feedback, either 
from star formation that is internal to the clouds 
themselves or from star formation activity in the 
nearby spiral arm. 

At smaller radii (i?<60”) there is relatively less 
massive star formation, and here the impact of 
shear and Coriolis forces due to disk differential 
rotation are better highlighted. These factors can 
lead to cloud dispersal and destruction at all radii 
in the inter-arm, but shear measured by the back¬ 
ground Oort A is notably larger here than beyond 
i?~70”, as shown in the top panel of Figure 
This measure of shear should be appropriate in the 
inter-arm, where we expect non-circular streaming 
motions to be negligible p] 

In the bottom panel ofFigure|^we also plot the 
radial profile of the ratio of the star formation rate 
surface density to the molecular gas surface den¬ 
sity in the inter-arm for reference. Beyond i?=60”, 
the high star formation per unit gas mass suggests 
that feedback may be more disruptive here than at 


smaller radii where, in contrast, shear is high. In 
what follows, we therefore distinguish between the 
“shear-dominated” and “feedback-dominated” re¬ 
gions of MSI’s inner disk, which we define as being 
galactocentric radii 41” <R<60” and 60” <i?<91” 
respectively. 

3.5. Phase changes and continuity in M51 

In Figures and we confirm that CO emis¬ 
sion supplies a compete picture of the evolution 
of MSl’s inter-arm molecular cloud population, 
i.e. that phase changes do not cause incomplete¬ 
ness in the PAWS catalog and the transition from 
molecular gas to atomic gas occurs over a much 
longer timescale than the orbital period (and the 
expected cloud lifetime). The left-most timescale 
in Fig. I^in particular represents the length of time 
over which the ISM is expected to remain molec¬ 
ular, calculated based on the requirement of con¬ 
tinuity between the atomic and molecular phases 
(i.e. the two dominant ISM components in M5I 
by mass). 

Follow ing iScoville fc Hersh|(|1979|) (and|Koda 
al. 20091, th 2 = thiMh^/Mhi where Mh^ a: 


et 

H 2 and 


"Mhi are the masses of molecular and atomic gas 
in the inter-arm (see Figure]^, ignoring the neg¬ 
ligible mass in ionized gas as well as conversion 
into stars (considering the long > 1 Gyr deple¬ 
tion timescale implied by the current rate of sty 


formation, ^2 AIq yr Schuster et al. 2009). 


The mass of molecular hydrogen m the mter-arm, 
shown in the top panel of Figure is calculated 
from the PAWS CO emission assuming a CO-to- 
H 2 conversion factor of A=2xl0^° cm“^ (K km 
s“^)“^, which has b een found to apply ac ross the 
PAWS field-of-view (Colombo et al. 2014a and ref¬ 
erences therein). The HI mass (also shown in Fig- 
ure[^ is estimated using the THINGS survey data 
(Wmer et al.|2008) assuming that the HI emission 
is optically thin . 

According to |Scoville & HershI (|1979|), the HI 
lifetime thi can be approximated by the shorter 
of either the dynamical time i'-^2ttravei here; see § 
^ or the spiral arm passage time tgp = 2Trfm{V,— 
Tip). We adopt the former as it is shorter than tgp 
between the m=2 spiral arms in M51, assuming 
the spiral pattern speed Up estimated by Quere- 
jeta et al. 2015 (in prep; see also Meidt et al. 
2013)(3 We obtain a more conservative estimate 
for th 2 by letting THi=ttravei as it everywhere 
underestimates the cycling timescale in the case 
where THi=tdyn ~ 2ttravei- This shorter timescale 
also allows for direct comparison between the esti¬ 
mated th 2 and the time window Uravei probed by 
the azimuthal span of our analysis region in the 
present inter-arm. 

The fact that the left bar depicted in Figure 
greatly exceeds Uravei therefore immediately sug- 


Note that in the spiral arms where streaming motions 
are larger, shear described by Oort A accounting for these 
motions (and including the background) can behave very 
differently, suggesting that the impact of shear on cloud 
stability may differ between the inter-arm and the arm. 


The spiral arm passage time would apply in the case 
that spiral arm passa ge prompts phase changes from atomic 
to molecula r gas (i.e. |Dobbs et al.|2008||Vazquez-Semadeni| 
let al .l2007ll and this happens taster than a dynamical time 
(i.e. [scdville fc Wilson|2004[l 
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gests that the disappearance of clouds from the 
PAWS catalog in the inter-arm is not the result of 
phase changes from molecular to atomic, and can 
instead be associated with a genuine finite cloud 
lifetime. As noted in § |2.3.2[ other galaxies may 
have lower molecular gas fractions, in which case 
the implied length of time that the gas stays in 
molecular form is reduced. Thus, we recommend 
first estimating the timescales described above be¬ 
fore applying our method in order to properly as¬ 
semble a multi-phase ce nsus o f clouds, where nec¬ 
essary. As suggested in § |2.3.2[ we also recommend 
assessing mass continuity across the inter-arm as 
demonstrated in the bottom panel of Figure in 
order to insure that these timescales are meaning¬ 
ful. 

In M51, H 2 traced by CO dominates the to¬ 
tal gas mass, which is almost equally distributed 
between zone I and zone II at all radii in the 
inter-arm. This suggests that the molecular phase 
largely captures mass continuity. The large varia¬ 
tions in the mass in cloud form from one side of the 
inter-arm to the other shown in the bottom panel 
of Figure therefore imply that clouds leave the 
population as a result of destruction rather than 
phase changes. 

Note that the above continuity argument also 
suggests that inter-arm clouds are converted di¬ 
rectly back to the atomic phase only very slowly, 
with a timescale which is again long compared to 
ttravei (third bar from the left), assuming once 
again THi=ttravei- This long timescale for con¬ 
version to atomic hydrogen implies that molecu¬ 
lar clouds in M51 would primarily evolve quickly 
back into their molecular hydrogen surroundings. 

The same continuity argument indeed implies a 
very fast conversion between clouds and the sur¬ 
rounding molecular gas (right bar). As noted ear¬ 
lier in § |2.3[ the cloud component of the inter¬ 
arm molecular gas in M51 (comprised of objects 
more massive than 10® M©) represents ^ 40% of 
the total inter-arm CO flux mapped by PAWS 
(see also Figure]^. This once again assumes that 
ttravei is the characteristic time for the molecu¬ 
lar (and atomic) gas to remain ‘outside’ clouds 
(i.e. assuming the medium is converted into clouds 
during passage through the spiral arm), and that 
the same CO-to-H 2 conversion factor applies to 
cloud and non-cloud CO em ission traced by PAWS 
(Liszt, Pety & Lucas 2010). Such a short cloud 
lifetime by this measure is consistent with the life¬ 
times found here (presented in § ffl, which also 
captures cloud conversion back to tne (molecular) 
non-cloud medium. But our estimate requires no 
prior knowledge of the timescale characteristic of 
the non-cloud medium, which may deviate from 
the assumed ttravei depending on the true cloud 
formation timescale. Our method is also indepen¬ 
dent of mass estimation, and so is free of the re¬ 


lated uncertainties, as discussed in § 2.4 


4. RESULTS 

4.1. Measurements and uncertainties 



Fig. 5. — Histogram of the average characteristic 
timescales for molecular hydrogen in the inter-arm of M51 
from throughout the zone 41”<i?<91”. The light gray bar 
shows the time for the ISM to remain molecular given con¬ 
version to atomic hydrogen, which is assumed to have char¬ 
acteristic lifetime rjji set by half the dynamical timescale 
{ttravei 1 fhe travel time from one of two spiral arms to the 
next; shown here as the white bar). The set of dark gray 
bars show the characteristic timescale for the part of the 
molecular hydrogen in cloud form, as cataloged by PAWS, 
assuming conversion either directly to atomic hydrogen or 
back to the non-cloud medium, including molecular hydro¬ 
gen. Radial profiles of the latter two timescales are shown 
in Figure]^ 


To best reveal the radial dependence of cloud 
dispersal and destruction processes expected in 
M51, we divide the inter-arm environment into a 
series of radial bins with uniform width. Each 
radial bin is further split at the midpoint of its 
azimuthal extent into two areas, zone I and zone 
II, in which we count the number of clouds Ni and 
Nil. Radial bins are discontinued at i?=90”=3.3 
kpc, the last radius in the field-of-view where it 
is possible to assign equal areas to zones I and II 
and cloud identification is reliable. The bin width 
was chosen to ensure sufficient numbers of clouds 
in each zone of each bin (a minimum of 4 and as 
many as 23) while at the same time matching our 
conservative estimate for the size of the maximum 
radial excursion expected for clouds. 

The uncertainties on our measurements reflect 
the changes to Ni and Nu when the boundaries 
of the inter-arm zones are modified. Part of the 
uncertainty is estimated by changing the location 
of radial bins by ±1” (the PAWS resolution; 15% 
of the bin width). This uncertainty is added in 
quadrature with the error arising from displace¬ 
ment of the inter-arm midpoint by 8°, accounting 
for the ^ 4 Myr uncertainty in our estimate of the 
inter-arm travel time. 

For reference, with this standard error budget 
the uncertainty on the total gas ma ss measured 


in either zone I or zone II is 30% (Colombo et 


al. 2014aI, comparable to the fraction of mass in 


cloud form. The significance of the mass-based 
diagnostics that include this error, which we use 
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Fig. 6. — (Top) Radial profiles of the mass in molec¬ 
ular gas (gray), atomic gas (gray dashed) and molecular 
clouds (black) in the inter-arm of M51. Meas urem ents 
and error bars adopt the scheme used later in |4.1| i.e. 
radially binned averages for which the uncertainties arise 
with modifications to the bin definition. Errors on the 
mass in clouds additionally include (and are dominated 
by) propagated measurement uncertainties on individual 
cloud masses. (Bottom) Ratio of the total (molecular plus 
atomic) gas mass in zone I compared to zone II (gray) and 
the ratio of the total mass in clouds in zone I compared to 
zone II (black). 


later in § |5.1| to interpret our measurements of 
Float and r, is therefore limited. However, the 
significance of our measurements of Float and r 
are not compromised. 

4.2. Evolution in M51s inter-arm cloud 
population 

In Figure we present our measurements of 
Float (right) and the cloud lifetime implied by this 
evolution in cloud numbers (left) for the inter-arm 
region of M51 probed by PAWS. Already by mid¬ 
way through the inter-arm almost 80% of the pop¬ 
ulation has been destroyed, implying a very short 
lifetime. Inside R ~ 70”, our method yields a 
characteristic cloud lifetime of only 20 to 30 Myr. 
Furthermore, we find only modest variation in the 
evolution of clouds numbers from zone I to zone 
II with galactocentric radius. This is in contrast 
to the increasing losses expected in a population 
with a uniform cloud lifetime as the travel time 
lengthens with radius. Instead, the trend on the 
right is qualitatively consistent with cloud disper¬ 
sal primarily through shear, which weakens with 
galactocentric radius and thus leads to longer life¬ 
times at larger radii (in the absence of other cloud 
destruction mechanisms). 

Quantitatively, moreover, we find that our life¬ 
time estimate agrees very well with the shear 
timescale (especially at R <70”; left), strongly 


suggesting that shear is the primary mechanism 
responsible for the finite lifetimes of clouds in 
M51’s inner disk. While the data are consistent 
with the shear trend over the range 40”<i?<85” 
within la, we can reject the constant average value 
of 31 Myr with 3a confidence. Note that even at 
the smallest radii where the inter-arm travel time 
is shortest, few clouds appear to survive in the 
presence of such strong shear. 

At larger radii (i? > 60”), though. Figure 
shows that the good agreement between the cloua 
lifetime and the shear timescale breaks down. At 
these larger radii, shear is slightly weaker and 
maps of star formation rate tracers (e.g. Ha, 
infrared, UV) suggest that star formation (and 
feedback) may have a greater influence on cloud 
evolution. Between 70”<i?<85”, the cloud num¬ 
ber density once again decreases from zone I to 
zone H, but not as many clouds are lost as would 
be expected under the influence of strong shear 
(i.e. Float falls below the shear prediction; right 
panel Figure [^. The cloud lifetime estimated by 
our method tnus appears lengthened in compari¬ 
son with the shear timescale (if only marginally). 
These trends are consistent with growth in the 
population between zone I and zone H at the 
same time as clouds are destroyed. Our measure¬ 
ments may therefore suggest that cloud forma¬ 
tion and/or transformation become increasingly 
important beyond R > 60”. As discussed in pj 
our method provides only an upper limit on the 
cloud lifetime unless the number of clouds added 
to the population is negligible compared to the 
number of clouds that are destroyed. Hence, our 
conversion from Float in the right panel of Fig¬ 
ure to the cloud lifetime shown in the left panel 
becomes highly uncertain at larger radii. In the 
next section, we examine the fate of clouds within 
M51’s inner disk in more detail, examining how 
their mode of destruction (shear or feedback) af¬ 
fects their fate (i.e. fragmentation/dispersal ver¬ 
sus transformation) and the lifetime that we esti¬ 
mate using our method. 

5. INTERPRETATION 
5.1. The fate of clouds 

Figure 1^ shows two additional diagnostic mea¬ 
sures of cloud evolution. In the top panel we show 
variation in the mass in clouds (black) and the 
mass outside clouds (gray) from zone I to zone H, 
calculated as {Mj — Mjj)/Mj. In the former case, 
the cloud mass in either zone i is measured as the 
sum of the masses of Ni individual clouds. In the 
latter case, the non-cloud mass in either zone is 
measured by subtracting the cloud mass from the 
total mass in gas in that zone, i.e. 

^non-cloud ~ ^gaedA^ — Mi ( 6 ) 

where At is the area covered by zone i and Egas in¬ 
cludes both atomic and molecular hydrogen traced 
by THINGS and PAWS. Note that cloud material 
returned to the non-cloud medium between zone 
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Fig. 7. — (Left) Radial profiles of cloud lifetimes measured for the inter-arm cloud population in M51. The solid black 
line shows our estimate using eq. [^and the measurements of Fi^st shown on the right. Error bars are propagated from the 
uncertainties on the measured Fiast- The black dotted (dashed) line shows the shear timescale A~^ (the inter-arm travel 
tim e ttravel)- The set of thick dashed gray lines show two independent cloud lifetime estimates based on continuity (see § 
\2A\ , with average values indicated by the two rightmost bars in Figure]^. (Right) Radial profile of the fraction of ‘lost’ 
clouds between zone I and zone II Fi^st- Error bars represent the change m lost fraction at each radius due to our standard 
uncertainty (see text for details). The dotted line represents the radial behavior in Fi^gt expected if the cloud lifetime is 
set by the shear timescale A~^, while the dashed line assumes a fixed lifetime of 40 Myr. Black arrows indicate the radial 
bin in which the number of clouds increase from downstream to upstream (see § |5.3|l . 


I and zone II is expected to rejoin the molecu¬ 
lar phase, as the timescale for conversion directly 
back to atomic hydrogen is much longer (see § 


2.3.2). 

In the bottom panel of Figure we show the 
change in the cloud-to-total mass ratio Rdouds be¬ 
tween zone I and zone II. Here, Rdouds in a given 
zone is the total mass in clouds divided by the 
total neutral (molecular plus atomic) mass in gas 
within that area. 

Measured quantities in each radial bin are shown 
with error bars that include the propagated uncer¬ 
tainties on cloud masses tabul ated in the PAWS 
catalog (Colombo et al. 2014a) and uncertainties 
aris ing from our standard error budget described 
in § 4.1[ While these errors can be quite large. Fig¬ 
ure ^ still provides a useful basis for interpreting 
the trends in Figure The two panels together 
suggest that clouds can undergo very different evo¬ 
lution depending on their radial location in the 
disk. 


5.2. Shear-dominated evolution 

To start, at small radii where shear is strong, 
there is signihcantly less mass in clouds in zone 
II compared to zone I suggesting that shear acts 
to disperse clouds (or cause a reduction in mass 
below our sensitivity limit). The bottom panel 
of Figure showing the evolution in the cloud- 
to-total mass ratio Rdouds from zone I to zone 
II tends to confirm this scenario. At ii<60”, 
the mean Rdouds/Rdouds is consistently (if only 
marginally) larger than one, suggesting that more 
of the total gas mass is in cloud form in zone I 
than in zone II. The inverse - more mass exterior 
to clouds (IO^Mq and higher) in zone II than zone 
I - is consistent with the idea that shear promotes 
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Fig. 8. — (Top) Radial profiles of the fractional mass in 
clouds lost from zone I to zone II in the inter-arm of M51 
(black) and the fractional mass in the non-cloud medium 
lost from zone I to zone II (gray). (Bottom) Radial profile 
of the ratio of the cloud-to-total mass in zone I R^, , 

clouds 

compared to the cloud-to-total mass in zone II (see 

text for details). 


mass loss and that this mass is returned to clouds’ 
surroundings. 

We find a hint of such an exchange between 
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clouds and the non-cloud medium in the evolu¬ 
tion in non-cloud mass from zone I to zone II, 
which tends to increase across the inter-arm at 
small radii (leading the gray line to drop below 
zero; top panel). We caution that since the cloud 
contribution to the total gas mass in each zone is 
on the order of its uncertainty, this is only weakly 
revealed by our measurements. 

5 . 3 . Feedback-dominated evolution 

The fate of clouds appears to be different at 
larger radii where shear is smaller and feedback 
has a potentially greater impact on clouds. Be¬ 
tween 60”<i?<85”, the loss of mass in cloud form 
between zone I and zone II is less than for re¬ 
gions at smaller galactocentric radii (Figure!^ top 
panel). Here as well the fraction of mass in aoud 
form (bottom panel) remains roughly fixed from 
zone I to zone II, as does the number of clouds. 
This suggests a scenario in which cloud destruc¬ 
tion is accompanied by the creation of new clouds 
from within the existing population, keeping the 
number of ‘lost’ clouds low. 

While a small Fiost might indicate little evolu¬ 
tion and a genuinely long cloud lifetime, this seems 
less likely given that the outer zone is very clearly 
impacted by star formation feedback (see Figure 
[^, which might be expected to limit cloud life¬ 
times below the inter-arm travel time. Instead, if 
the population contains growth at a rate similar 
to the rate of destruction, the number of clouds 
would remain fixed and eq. would overestimate 
the cloud lifetime. Below we consider several sce¬ 
narios in which cloud destruction might be bal¬ 
anced by cloud formation. 

5 . 3 . 1 . Transformation via feedback 

In the simplest case, feedback might act to split 
or merge clouds, transforming the cloud popu¬ 
lation, rather than completely destroying or dis¬ 
persing clouds. New clouds would thus emerge 
from within the existing population, as descen¬ 
dants of transformed clouds. Note that if the cloud 
population is undergoing transformation. Figure 
[^suggest that clouds must be merging as well as 
splitting, as splitting alone would be expected to 
increase the number of clouds by zone II. How¬ 
ever, we cannot rule out that some clouds may 
also be completely destroyed, e.g. in the man¬ 
ner described in the previous section; mass loss 
might follow from shear, which is non-zero at these 
radii, or through the process of star formation it¬ 
self, which might consume a large fraction of the 
cloud mass. In some models, the star formation 
efficiency per free-fall time can be as high as 0.2- 


0.3 ( 

e.g. Klessen & Burkert 2000 Bonnell et al. 

2003 

IVazquez-Semadeni et al. 20031 Clark et al. 

2005 

1 . 


5 . 3 . 2 . Population growth 

It is also possible (if perhaps less likely) that 
new clouds at these radii are formed completely 
independently of the destruction process. (In fact. 


we allow for this possibility at all radii, with the 
caveat that our measure of r is an upper limit 
on the true lifetime.) To account for the trends 
between 60”<i?<85” in Figure cloud creation 
would have to roughly balance destruction and 
the genuine cloud formation timescale Tg^ow in the 
population would need to be comparable to the 
cloud lifetime. 

Although genuine population growth in the 
inter-arm may be surprising, cloud formation may 
not always be limited to the spiral a rms. In the 
simulations of |Dobbs fc Pringle |2013[ spirals that 
radially decrease in strength become less impor¬ 
tant as a site of cloud formation. In M51, the 
strong two-armed spiral transit ions to a weaker, 
material-like pattern a t i?=85” (Meidt et al.|2013 
Colombo et al.|[^14b ) 

i'he importance of inter-arm cloud formation 
might moreover be expected to increase with ra¬ 
dius, since the time between arms is longer at 
larger radii. This could provide sufficient time to 
form new clouds from the material that was re¬ 
turned to the ISM by previously destroyed inter¬ 
arm clouds. Clouds could also form directly from 
the non-cloud medium, presumably at a fixed rate, 
thus leading our measurements to overestimate 
the true cloud lifetime progressively more with ra¬ 
dius. To confirm whether clouds can form in the 
inter-arm, larger, more spatially extended molec¬ 
ular cloud surveys are necessary. At present, we 
take our measurements in the zone 60” <i?<83” as 
most likely representing transformation within the 
population stimulated by feedback, as discussed in 
the previous section, rather than cloud destruction 
accompanied by independent population growth. 

5 . 4 . (In)sensitivity to level of virialization, mass 
and surface density 

The strength (or weakness) of variations in the 
cloud population from zone I to zone H at a given 
radius could depend on whether all clouds expe¬ 
rience the same limit to their lifetimes, or are 
equally susceptible to destruction. The stability 
of clouds against dispersal or destruction could de¬ 
pend on, e.g., cloud mass or the balance between 
internal kinetic energy and gravitational potential 
energy (level of virialization). We might therefore 
expect that some clouds are never dispersed while 
others may evolve more rapidly than the rest of a 
given population. 

In MSI, however, we find no significant change 
in average clo ud properties, including th e virial 
parameter a (Bertoldi & McKee 19921, cloud 
mass or surface density, from zone I to zone H 
and the two pop ulations as a whole app ear very 
similar (see also Colombo et al. 2014a). Addi¬ 
tional tests moreover show no strong link between 
cloud lifetime and any given property. Specifically, 
we separated clouds into subpopulations and then 
compared the number of clouds in zone I with a 
given property to the number of clouds in zone H 
with that same property at each radius. The mea¬ 
sured Fiost and lifetime r in each subpopulation 
show little significant difference, e.g. the evolu- 
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tion in the number of clouds with a<2 is nearly 
indistinguishable from that of the subpopulation 
with a>2. While it thus appears that no sub¬ 
set of clouds is more susceptible to a particular 
destruction mechanism than any other, we em¬ 
phasize that additional splitting of the inter-arm 
populations in zone I and zone II likely leaves in¬ 
sufficient cloud statistics to reliably apply eq. ([^. 

6. DISCUSSION 

6.1. Cloud destruction processes 

In the previous section we used the combina¬ 
tion of relatively little massive star formation and 
strong shear at i?<60”, and the presence of more 
massive star formation at radii i?>60” where shear 
is weaker, to isolate the influences of shear and 
star formation feedback on clouds. We emphasize, 
however, that we do not directly observe cloud 
destruction, and thus cannot absolutely conclude 
which process or processes limit cloud lifetimes. 
So although we do find radial variation in the way 
the cloud population evolves across the inter-arm 
- which we interpret as the signature of two main 
modes of cloud destruction - we can attribute this 
only generally to shear and feedback. 

Indeed, in the present study we have very little 
leverage on the way the two processes may work 
together. The susceptibility of clouds to shear 
could still depend on the impact of star formation, 
which we observe to be non-zero even at i?<60” in 
the inter-arm (see Figure]^. Star formation feed¬ 
back could be necessary to first destabilize clouds, 
or it could serve to enhance the dispersal process. 
Likewise, even when shear is relatively weak, it 
may still assist in cloud destruction through other 
processes, including feedback. (Note, though, that 
the fate of clouds at large radii seems inconsis¬ 
tent with evolution driven entirely by shear, as 
the overall mass in clouds does not decrease from 
zone I to zone II; see Figure]^. 

6.2. Cloud lifetime measurements 

Figure reveals a close relation between the 
shear timescale measured by A~^ and the cloud 
lifetime at i?<60”. This suggests that, even if star 
formation (and its associated feedback) is required 
as a prior source of cloud destabilization, when 
shear is strong enough it takes over and sets a more 
important limit on the cloud lifetime r, indepen¬ 
dent of the strength or pattern of feedback. The 
lifetimes of all clouds (at least those more massive 
than 10 ®Mq as considered here) seem equally lim¬ 
ited by the shear timescale when shear is strong. 

ED 

But as shear weakens and the shear timescale 
increases with galactocentric radius, at some lo¬ 
cation in the disk we expect destruction via feed¬ 
back to take over and set the cloud lifetime. Ac- 

In the less likely alternative scenario (introduced in 
1 ]5.3[| , in which cloud destruction is accompanied by inde¬ 
pendently forming new clouds at all radii, our measure¬ 
ments could imply even shorter cloud lifetimes than shown 
in Figure]^ 


cording to Figure we would argue that feed¬ 
back is accompanira by transformation (creation 
and destruction). Our measurements of t in Fig¬ 
ure [ 7 ] thus provide only an upper bound on the 
cloud lifetime due to feedback and not a direct es¬ 
timate. Yet we can extract an estimate for the 
cloud lifetime due to feedback using the fact that 
the associated transformation keeps -Fiost in the 
population low. The location where this signature 
emerges marks the location where the lifetime due 
to feedback is shorter than the shear timescale. 
We can thus estimate the feedback timescale from 
the shear timescale at this radius. 

In M51, a sustained drop in Fiost occurs at 
i?^70” where 30Myr, suggesting that the 

lifetime due to feedback is around rwdOMyr. 
(This location is not surprisingly very close to the 
radius R—60” that, by construction, distinguishes 
between zones dominated by either shear or feed¬ 
back.) Eqs. ^ and ([^ can provide an alter¬ 
native estimate, assuming a particular model for 
how the transformation via feedback occurs (e.g. 
equal-mass splitting) to set the relation between 
the rates of cloud creation and destruction. 


6 . 3 . Implications of short CMC lifetimes 
6 . 3 . 1 . Cloud evolution and star formation 


The cloud lifetimes measured with our technique 
agree very well with the few existing observation- 
based estimates (typically t«20-30 Myr) made 
with a completely independent method (i.e. link¬ 
ing c louds with stellar clusters at various young 
ages; [Bash et al.|1977[|Leisawitz et al.|I989[[Miura 

et al.j 2012 Kawamura et al. 2009|). Such short 

lifetinies agree with the picture of rapid cloud evo¬ 
lution that emerges from numerical simulations, 
where clou ds typically have short (10-20 Myr) life¬ 
times (e.g. Dobbs et al. 2011 Dobbs et al. 


Dobbs & Fringle 2013 b We can thus conhrm 

that, whether in star-forming disk galaxies or in 
low-mass systems, clouds have sufficiently short 
lifetimes that they are disrupted after a few free- 
falls times, as previously suggested (Elmegreen 


2000[ Ballesteros-Paredes fc Hartmann|2U07[ |Mur- 

ray||2011 ). 


Th this light, cloud longevity would appear to 
provide an unsatisfactory resolution to the is¬ 
sue of low observed star form ation efficiencies 
(e.g. Krumholz fc Tan] 2007). Instead, short¬ 


lived clouds may suggest that only a fraction of 
the cloud population undergoes ‘active star for¬ 
mation’. Clouds can then be individually effi- 


els; IKlessen & BurkertI 

20001 Bonnell et al. 1120031 

Vazquez-Semadeni et a 

r 2003 Clark et al. 1200^ 

even while the overall e 

thciency in the current gas 


reservoir remains low. 


6 . 3 . 2 . The cycling of molecular gas from diffuse to 
bound objects 

Our measurements further test the notion of 
cloud longevity by directly contradicting several 
arguments invoked in its favor. We find that 
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clouds more massive than 10^ Mq in M51 sur¬ 
vive for much less time than the characteristic 
timescale of the molecular phase (see Figure 
previously upheld as a measure of the cloud life¬ 
time. The difference in timescales immediately 
suggests that a significant fraction of the molec¬ 
ular gas exists outside of clouds with masses 10^ 
Mq or more. Indeed, roughly half of PAWS CO 
emission has been determined to be in the form 
of an extended component, rather than in com¬ 


pact s tructures (Pety et al. 2013 Colombo et al. 


2014a). 

T he existence of a molecular non-cloud medium 
itself suggests that gas can stay molecular all 
the while clouds are dispersed and reassemble. 
Comparing the mass in clouds with the mass in 
M51s non-cloud medium, we find that the con¬ 
version of the cloud component can be fast. We 
independently estimate a very short characteris¬ 
tic cloud t imesca le based on the continuity argu¬ 
ments in § 2.3.2 i.e. tcont—^clouds/^^non — cloud ^ 
Tnon-cioud where Mclouds is the mass in cloud 
form, Mnon-cioud IS the (molecular and atomic) 
mass outside clouds and Tnon-cioud is set to travel ■ 
This timescale, which is plotted in Figure [7j is 
comparable to the cloud lifetimes measured with 
our framework. (The ave rag e is represented by the 
rightmost bar in Figure The difference may 
indicate that our measurements overestimate the 
true lifetime. But it more likely signifies that gas 
remains outside clouds, in t he non -cloud medium, 
longer than we adopted in § 2.3.2[ 

Recall that we conservatively let 
'^non—cloud — ^travei'^^orb/^^: likely Underestimating 
the actual time it takes for the non-cloud medium 
to be converted into clouds by a factor of 2. If we 
take our measurements as the true lifetime, then 
comparing r with tcont suggests that Tnon-doud 
may be closer to 2r « 2A~^ < 2ttravei, i-e. the 
timescale for cloud formation from the non-cloud 
medium is on the order of a dynamical time. 

If the reverse is true, and tcont is a more real¬ 
istic measure of the cloud lifetime than our esti¬ 
mate, then eq.([^ implies that the inter-arm cloud 
population contains growth characterized by a 
timescale that is roughly 2tcont=‘^ttravei- 

Interestingly, in either case, to match tcont and 
r in Figure |r requires that clouds leave the inter¬ 
arm cloud population quicker than they can be 
replenished from the non-cloud medium. Conse¬ 
quently, the cloud-to-non-cloud mass ratio should 
vary as a function of azimuth, from one side of the 
galaxy to another. As we do not observe this vari¬ 
ation, we can conclude that, as expected, clouds 
form in the spiral arm and these clouds feed the 
inter-arm population. To maintain the inter-arm 
cloud-to-non-cloud mass ratio over the course of 
one orbital period (from 0 to 27r), as observed, the 
formation timescale for clouds in the arm must be 
on the order of 2r. 

These considerations lead us to conclude that 
(large-scale) galaxy dynamics regulates the cy¬ 
cling of molecular material between a diffuse state 
and a more bound state (represented by ‘clouds’) 


from which stars can form. Clouds emerge from 
the diffuse state, and their recognizability might 
depend not only on the molecular content of a 
galaxy disk but also on its dynamical character, 
i.e. the total (baryonic) mass and its distribu¬ 
tion, including (non-axisymmetric) bar and spiral 
structures. 

6.4. Implications for other cloud populations 

Our cloud lifetime measurements in M51 not 
only confirm, but also extend, the picture of short 
cloud lifetimes so far supported by only few exist¬ 
ing measurements. By probing cloud populations 
under sufficiently different conditions than previ¬ 
ously considered (spanning a larger range in radius 
in the disk of a spiral galaxy) we establish a ba¬ 
sis for understanding how cloud lifetimes may be 
expected to vary in general. 

Here we found that everywhere the shear 
timescale is faster than the feedback timescale 
the inter-arm cloud lifetime is set by A“^, with 
very little sensitivity to the rate (or amount) 
of star formation. We therefore propose that, 
in other circumstances, knowledge of the shear 
timescale should be sufficient to predict the 
cloud lifetime. Whereas feedback arguably 
proceeds with a universal timescale set by the 
~ 30 Myr lifetimes of massive OB stars, the 
shear timescale, which can drop below 30 Myr, 
determines where feedback becomes the dominant 
limit to the cloud lifetime. Elsewhere, shear it¬ 
self can provide a direct constraint on the lifetime. 

In Spiral Arms 

In spiral arm environments, for example, where 
characteristic strong streaming motions can 
locally reduce shear, the increase in the shear 
timescale suggests that spiral arm clouds can 
be longer-lived than their counterparts in the 
inter-arm. 


In Late-type Galaxies 

Shear should be a similarly good predictor of 
cloud lifetimes throughout cloud populations, 
particularly in more massive disks. Since the 
shear timescale depends on the shape and maxi¬ 
mum of the circular velocity in the disk, it will 
likely remain shorter than the feedback timescale 
throughout most of all but the lowest mass 
galaxies. Cloud lifetimes would be expected to 
increase with decreasing galaxy mass, until shear 
timescale overall exceeds the feedback timescale. 
This arguably explains why clouds in the LMC 
have comparably sh ort lifetimes (20-30Myr; 


Kawamura et al. 2009) as clouds in M51, despite 
the much longer shear timescale. 


In Early-type Galaxies 

The high shear rate in the centers of massive, 
early-type galaxies may increase the likelihood 
that some molecular gas may never have the 
chance to form stars, if cloud lifetimes (set by the 
shear timescale) become comparable to, or shorter 
than, the free-fall time. This might lead to less 
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efficient star formation and lengthened gas deple¬ 
tion times, such as measured in early-type galax¬ 
ies with CO detections, where depletion times 
typically exceed those measured in normal star¬ 
formi ng galaxies by a factor of 2.5 (Davis et al.| 

Mil. 


This provides a more compelling interpretation 
for the lengthened r^ep in such galaxies than just 
the shape of the rotation curve (as suggested by 
Davis et al.||2014 |, considering that 




TT 1 — /3 


(7) 


where the rotation curve shape is parameterized 
as /3=d In V/d In R. While eq. 0 suggests 
that short A~^ may lead to longer Tdp.n in galax ies 
where /3 is low, as found by|Davis et al.| (|2014|), it 
is also clear that a low (3 does not always imply 
a s hort A~^. This reso lves the discrepancy noted 
by Davis et al. (20141 that the r^ep in the disks 
of nearby late-type galaxies, where /3~0 but tort 
is long, are not comparably long as in early type 
galaxies where torb is short, but /3^0 as well. 


7. CONCLUDING REMARKS 

In this paper we present a new technique to mea¬ 
sure the lifetimes of giant molecular clouds that 
tracks formation and destruction within cloud 
populations through cloud number statistics. Our 
framework uses the travel time between spiral 
arms as a hducial clock, rather than a star 
formation-related timescale, and yields a charac¬ 
teristic cloud lifetime estimate even when cloud 
masses are poorly determined (due to e.g. uncer¬ 
tainties in the CO-to-H 2 conversion factor). As 
these ambiguities are avoided, cloud lifetimes mea¬ 
sured with our technique can yield unique insight 
into the dynamical influences on clouds, using only 
the information that is readily accessable in cur¬ 
rent and future molecular cloud surveys. 

Our hrst application of the method leverages 
the large cloud population in the inter-arm of 
M51 surveyed by PAWS. In a series of radial bins, 
we relate the change in cloud numbers across the 
inter-arm to an estimate of the average cloud 
lifetime at that radius. We use our detailed 
knowledge of gas dynamics and massive star 
formation across the PAWS field-of-view to iden¬ 
tify two radial zones where cloud destruction is 
dominated by shear and star formation feedback. 


respectively. Our analysis suggests the following 
conclusions: 

1. GMC lifetimes in the inner disk of M51 are 
short, typically 20 to 30 Myr. 

2. Shear due to galactic differential rotation is the 
primary limit to cloud lifetimes when the shear 
timescale is shorter than the feedback timescale 
(^30 Myr, the average lifetime of OB stars). 

3. At galactocentric radii in M51 where we expect 
shear to be the dominant mode of cloud disrup¬ 
tion, the evolution in the mass and number of 
GMCs across the inter-arm indicates that clouds 
are effectively dispersed. At larger galactocentric 
radii where there are strong signatures of high 
mass star formation, the cloud population appears 
to undergo transformation (i.e. an exchange of 
mass between high- and low-mass clouds) rather 
than complete dispersal. 

4. Since shear depends on the shape and maxi¬ 
mum amplitude of the galaxy’s circular velocity 
curve, low mass disks should contain longer-lived 
cloud populations than higher mass disks. In 
the most massive systems and in the concen¬ 
trated centers of galaxies, where short shear 
timescales approach the free-fall time, clouds may 
be so short-lived that star formation is suppressed. 

5. Clouds in M51 are shorter-lived than the 
characteristic lifetime of molecular hydrogen, 
implying that molecular material can continually 
cycle between clouds and their (non-cloud) sur¬ 
roundings. Based on our results, we suggest that 
conversion from a diffuse molecular phase into 
bound objects is regulated by large-scale galaxy 
dynamics. 

Future applications of our method to high- 
resolution, wide-held CO surveys of galaxies with 
a diverse range of Hubble types will be essential for 
testing these conclusions and conhrming the role 
of short GMC lifetimes in regulating the star for¬ 
mation efficiency observed across a range of spatial 
scales. The appropriate datasets to test our model 
are available with the advent of regular science op¬ 
erations at ALMA. 
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